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摘 要 :， 现 有 UT1-UTC 预 报 模式 在 进行 周期 项 与 残 差 项 拟 合 分 离 时 ， 通 常 没 有 考虑 最 小 二 
乘 拟 合 序 列 的 端 部 畸变 现象 (数据 处 理 中 称 为 端 部 效应 ) ， UES RAMS. EXT 
最 小 二 乘 拟 合 存 在 的 端 部 畸变 现象 ， 首 先 采 用 时 序 分 析 方 法 在 UT1-UTC 序 列 两 端 渤 了 数据 延 
拓 ， 形 成 一 个 新 序列 ， 然 后 用 新 序列 求解 最 小 二 乘 外 推 模型 系数 ， 最 后 再 0 合 最 小 二 乘 外 推 
模型 及 神经 网 络 对 UT1-UTC 序 列 进行 预 则 。 结 果 表 明 ， 在 UT1-UTC 序 列 端 部 增加 延 拓 数据 ， 可 
以 有 效 地 抑制 最 小 二 乘 拟 合 序 列 的 端 部 畸变 ; 相对 于 常规 最 小 二 乘 外 推 模型 ， 基 于 端 部 效应 
改善 的 最 小 二 乘 (edge-effect corrected least squares, ECLS) 外 推 模型 的 UT1-UTC 中 
KH TH Ek SAA 
关键 词 : UT1-UTC 预 报 ; 预报 模型 ; 最 小 二 乘 ; 端 部 效应 ; 端 部 延 拓 
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地 球 自转 运动 可 以 用 地 球 定 向 参数 (Earth Orientation Parameters, EOP) 表征 。 地 球 定向 


参数 包括 岁差 、 章 动 、UT1-UTC 及 极 移 的 两 个 分 量 X, 和 7》,。 地 球 定向 参数 是 实现 地 球 参考 


> 


坐标 系 与 天 球 参考 坐标 系 相 互 转换 的 必需 参数 ， 在 天 文 地 球 动力 学 研究 、 深 空 探测 及 卫星 


航 等 领域 有 重要 应 用 。 甚 长 干涉 基线 测量 (Very Long Baseline Interferometry, VLBI). 4 %4) 


卫星 定 轨 和 无 线 电 定位 C Doppler Orbitography and Radiopositioning Integrated by 


Satellite, DORIS) 和 全 球 卫 星 导 航 系 统 (Global Navigation Satellite System, GNSS) 等 空间 
测 地 技术 是 测量 地 球 定向 参数 的 主要 手段 ， 其 中 UTI-UTC 和 极 移 的 测量 精度 可 以 分 别 达 到 
10ps 和 100hasi， 但 复杂 的 资料 处 理 过 程 致使 地 球 定向 参数 的 获取 存在 一 定 的 延迟 。 由 于 
深 空 探测 等 领域 对 地 球 定 向 参数 实时 测量 值 有 重要 需求 ， 而 卫星 自主 导航 又 对 地 球 定向 参 
数 中 长 期 预报 值 有 重要 需求 ， 因 此 对 地 球 定向 参数 进行 高 精度 的 短期 与 中 长 期 预报 是 非常 
必要 的 。 在 地 球 定向 参数 的 几 个 分 量 中 ，UTI1-UTC 是 变化 最 快 、 最 难 预报 的 一 个 分 量 。 

目前 有 多 种 UT1-UTC 预 报 方法 ， 其 中 多 数 方法 结合 最 小 二 乘 (least squares, LS) 外 推 和 
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其 它 模 型 进行 UT1-UTC 预 报 * ， 这 些 方 法 的 思路 为 ， 首 先 采 用 最 小 二 乘 外 推 模型 提取 
UTI-UTC 序 列 中 的 周期 性 成 分 ， 并 对 其 进行 外 推 ， 然 后 利用 神经 网 络 (Neural Network, 
NMN) 或 者 自 回归 (Auto Regressive, AR) 模型 等 对 最 小 二 乘 拟 合 残 差 进行 建 模 、 预 报 ， 最 后 再 
将 周期 项 和 残 差 项 的 外 推 值 相 加 得 到 UT1-UTC 预 报 值 。 在 实际 应 用 中 发 现 ， 在 利用 最 小 二 
乘 外 推 模型 对 UT1-UTC 观 测 资料 进 行 拟 合 时 ， 在 拟 合 序列 的 两 端 存在 发 散 畸 变现 象 ， 这 种 
现象 在 数据 处 理 中 称 为 端 部 效应 。 端 部 效应 使 残 差 项 与 周期 项 的 预报 值 出 现 偏 差 ， 最 终 


导致 UT1-UTC 预 报 值 不 准确 。 本 文 针 对 UT1-UTC 预 报 中 最 小 二 乘 拟 合 出 现 的 端 部 畸变 现象 ， 


在 采用 最 小 二 乘 外 推 模型 对 UT1-UTC 序 列 进行 拟 合 之 前 ， 先 利用 端 部 延 拓 方 法 对 UT1-UTC 
观测 资料 进行 数据 延 拓 ， 即 在 UT1-UTC 序 列 的 两 端 增加 应 用 时 序 分 析 方 法 延 拓 出 的 若干 数 
据点 ， 形 成 一 个 新 序列 ， 然 后 用 新 序列 求解 最 小 二 乘 外 推 模型 系数 ， 最 后 再 基于 最 小 二 乘 
外 推 模型 对 UT1-UTC 序 列 中 的 周期 项 进行 外 推 。 数 值 分 析 表 明 ， 通 过 在 UT1-UTC 观 测序 列 
端 部 增加 统计 延 拓 数据 ， 可 以 有 效 地 抑制 端 部 效应 ， 从 而 改进 UT1-UTC 预 报效 果 。 


1 预报 方法 


1.1 数据 预 处 理 
本 文 UT1-UTC 观 测 资料 来 源 于 国际 地 球 自 转 与 参考 系 服务 (International Earth Rotation 


ar 


and Reference Systems Service, IERS) 发 布 的 EOP-08-C04 序 列 ， 其 中 UT1-UTC 数 据 的 采样 


间隔 为 1d。UTI-UTC 观 测 资料 中 含有 闽 秒 (Leap Second) 及 多 种 周期 项 、 准 周期 项 ， 其 中 ， 
对 于 周期 为 5d 一 18. 6a 的 62 个 固体 地 球 带 谐 潮 汐 项 应 用 IERS 协 议 CERS Conventions) 给 出 
的 经 验 公式 予以 扣除 外， 扣除 62 个 固体 地 球 带 谐 潮 汐 项 后 的 UT1-UTC 称 作 UTIR-UTC， 然 
后 再 去 掉 UT1R-UTC 序 列 中 闽 秒 ， 获 得 UTI1R-TAI 序 列 。 本 文 对 UT1-UTC 的 预测 实质 上 是 针 


对 UTIR-TAI 的 预测 。 
1. 2 最 小 二 乘 外 推 模型 

UT1-UTC 序 列 扣 除 62 个 固体 地 球 带 谐 潮汐 项 后 ， 还 含有 长 期 趋势 项 、 周 年 项 、 半 周年 项 
等 周期 性 变化 成 分 ， 长 期 趋势 项 与 周年 项 、 半 周年 项 利用 如 下 模型 进行 拟 合 、 外 推 : 


2 
forra =a+bt+ —[¢,cos(2nt/T,) +d, sin(2zt /T,)] 


i=l 


Ko, D Tankaha PHA a EM ee A, WATI =182.62d, 、 =365.24d 


Ml a, DRREAMEM, c, d 代表 半 周 年 项 参数 ，c, ， d, 代表 周年 项 参数 。 这 些 未 


知 参数 可 以 根据 最 小 二 乘法 求解 。 
1.3 神经 网 络 模型 


反 向 传播 (Back Propagation, BP) 神经 网 络 是 当前 应 用 较为 广泛 的 神经 网 络 算法 之 一 ， 


但 这 种 算法 存在 收敛 速度 慢 、 过 拟 合 及 泛 化 能 力 弱 等 不 足 。 鉴于 此 ， 本 文选 用 近年 来 发 展 起 
来 的 极限 学 习 机 Cextreme learning machine, ELM) 神经 网 络 算法 对 最 小 二 乘 拟 合 残 差 序列 
进行 建 模 、 预 报 。 与 反 疝 传播 神经 网 络 算法 相 比 ， 极 限 学 习 机 算法 训练 速度 非常 快 ， 不 需 

人 为 设置 网 络 参数 ， 且 具有 很 强 的 泛 化 能 力 “， 极 限 学 习 机 神经 网 络 在 地 球 定向 参数 预报 
中 已 经 取得 成 功 应 用 "4 ， 其 原理 如 下 所 述 。 


给 定 N 个 训练 样本 (x,, y,), ， 其 中 输入 向 量 x =[x,.x,.-5x, ] R”， 输 出 向 量 


Vi =YV Yn] R”，m，n 分 别 表示 输入 向 量 维 数 与 输出 向 量 维 数 ， 则 具有 六 


个 隐 含 层 节 点 的 单 隐 层 前 馈 神 经 网 络 能 够 以 零 误差 逼近 这 N 个 样本 ， 即 有 下 式 成 立 : 


h . 
yp X 8( 19; b,), i=1,2,---,N 


(2) 


Hh, w,=[Wy.Wasts Win] 表示 第 j 个 隐 含 层 节 点 和 输入 节点 之 间 的 连接 权 值 ， 


B, =[Bn,Bjs，…, Bn ] 表示 第 j 个 隐 含 层 节点 和 输出 节点 之 间 的 连接 权 值 ，b) 表示 第 J 


个 隐 含 层 节点 的 偏 置 ，g(X) 表示 激活 函数 。(2) 式 可 以 表示 为 矩阵 相 乘 的 形式 : 


Y=Hp, 


(3) 
RF, HARRA hE RE: 


ame, 的 ) = g(w, 六) 


a 


H= 


alr 


(w, Gey 的 ) a gw, Xy b,) Nh 
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(4) 


PoE 
B= J Y= F (5) 


极限 学 习 机 神经 网 络 的 训练 步骤 如 下 : CG) BEALE A BU Ww). E Db, 


j=1,2,…, 有 hh; (2) +RRA RHE A; (3) 计算 输出 层 权 值 p 古 了 i ， 其 中 


Hi 为 矩阵 H 的 Moore-Penrose 广 义 i. 

> 1.4 预报 过 程 
| 由 于 最 小 二 乘 外 推 和 神经 网 络 组 合 模型 (LS+NN) 不 仅 具 有 较 高 的 预报 精度 ， 还 具有 
= 较 好 的 长 期 预测 稳定 性 ， 故 基于 LS+NN 模型 的 UTI-UTC 短期 预测 值 精度 较 高 ， 可 以 很 好 
一 地 反映 UTI-UTC 的 变化 规律 ,作者 曾 将 端 部 数据 延 拓 的 思想 和 方法 应 用 于 日 长 变化 资料 
=D 分 析 中 ， 在 改善 最 小 二 乘 拟 合 的 端 部 效应 方面 进行 了 研究 ， 取 得 了 较为 明显 的 效果 “0。 本 
m 文 尝试 将 端 部 数据 延 拓 的 思想 和 方法 应 用 于 UT1-UTC 预报 ， 以 限制 UTI-UTC 资料 拟 合 中 
最 小 二 乘 拟 合 出 现 的 端 部 畸变 现象 。 
过 考虑 最 小 二 乘 拟 合 端 部 效应 的 UTI-UTC 预测 方法 和 常规 方法 的 差别 之 处 在 于 ， 该 方法 
fi. 在 对 UTI-UTC 观测 序列 建立 趋势 项 及 周期 项 最 小 二 乘 外 推 模型 之 前 ， 先 应 用 统计 学 方法 在 
UTIR-TAI 序列 首部 和 尾部 进行 数据 延 拓 ， 以 抑制 最 小 二 乘 拟 合 端 部 畸变 ， 预 报 过程 如 下 : 

(1) 首先 通过 (1) 式 对 UTIR-TAI 序 列 作 最 小 二 乘 拟 合 ， 建 立 周期 项 及 趋势 项 外 推 
模型 ， 然 后 利用 神经 网 络 模型 对 最 小 二 乘 拟 合 残 差 序列 进行 建 模 、 预报， 最 后 联合 最 小 二 乘 
外 推 模型 和 神经 网 络 模型 (LS+HNN) 在 UTIR-TAI 序 列 首部 和 尾部 分 别 外 推荐 干 数据 点 ， 
这 样 UTIR-TAI 序 列 加 上 首部 和 尾部 外 推 的 若干 数据 点 就 构成 了 一 个 新 序列 。 

(2) 利用 新 序列 求解 最 小 二 乘 外 推 模型 系数 ， 即 用 新 序列 重新 建立 趋势 项 及 周期 项 最 
小 二 乘 外 推 模型 ， 然 后 再 结合 最 小 二 乘 外 推 模 型 和 神经 网 络 模型 对 UTIR-TAI 序 列 作 预 报 。 

(3) 将 周期 为 54 一 18. 6a AY 62 个 固体 地 球 带 谐 潮 汐 项 及 头 秒 恢复 到 UTIR-TALI 序列 
的 预测 值 中 即 可 得 到 UT1-UTC 预测 值 。 


2 数值 分 析 


a 


chinaXiv:201805.00153v1 
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选取 2001. 1. 1 一 2016. 9. 20 期 间 的 UT1-UTC 观测 序列 进行 数值 分 析 ， 其 中 2010.1.1 
一 2016. 9. 20 为 预测 期 ， 建 模 数 据 长 度 为 10 年， 每 隅 7d 预 报 1 次 ， 总 共 进 行 了 300 次 预测 。 
图 1(a) 、(b) 分 别 给 出 了 2001. 1. 1 一 2016. 9. 20 时 段 的 UT1-UTC 测量 序列 及 其 扣除 闽 秒 和 


62 个 固体 地 球 带 谐 潮汐 项 的 UT1R-TAI 序 列 。 


为 了 验证 端 部 延 拓 方 法 对 最 小 二 乘 拟 合 端 部 效应 的 改善 


效果 ， 首 先 比较 端 部 延 拓 前 后 


最 小 二 乘 拟 合 效果 。 图 2 绘 出 了 2001 年 1 月 1 日 至 2010 年 12 月 31 日 时 段 UTIR-TAI 序 列 
及 其 端 部 数据 延 拓 点 ， 其 中 ， 端 部 延 拓 数 据点 数 为 360， 即 


点 。 从 图 


12 月 31 


在 首尾 两 端 各 延 拓 180 个 数据 


2 可 以 看 到 ，UTIR-TAI 端 部 延 拓 数据 点 与 原始 序列 符合 的 较 好 ， 可 以 很 好 地 反映 
原始 UTIR-TAI 序列 的 变化 规律 。 图 3 绘 出 了 端 部 数据 延 拓 前 后 2001 年 1 月 1 日 至 2010 年 
时 段 UTIR-TAL 序列 的 最 小 二 乘 拟 合 序列 ， 基 于 端 部 效应 改善 的 最 小 二 乘 (Edge- 


effect Corrected Least Squares, ECLS) 表示 端 部 效应 修正 的 最 小 二 乘 拟 合 序列 。 为 了 更 加 直 


数据 延 拓 方 法 对 最 小 二 乘 拟 合 端 部 畸变 的 改善 效果 ， 对 图 3 作 局 部 放大 


Ca 


合 残 差 数据 点 ， 图 4(b) 22H 


岗 、 清 楚 地 展示 端 部 
bE, Fl 4(a) 绘 出 了 端 部 数据 延 拓 前 后 最 小 二 乘 拟 合 序 列 首部 前 180 个 历 元 的 最 小 二 乘 拟 
上 了 端 部 数据 延 拓 前 后 最 小 二 乘 拟 合 序列 尾部 最 后 180 个 历 元 的 


最 小 二 乘 拟 合 残 差 数据 点 。 


从 图 


3 和 图 4 可 以 发 现 ， 与 直接 对 原始 UT1IR-TAI 序 列 进行 拟 合 相 比 ， 端 部 数据 延 拓 


后 最 小 二 乘 拟 合 残 差 序 列 在 首部 和 尾部 更 加 接近 于 0， 换言之， 端点 延 拓 后 最 小 二 乘 拟 合 
的 UTIR-TAI 序 列 在 首部 、 尾 部 和 原始 序列 吻合 的 更 好 ， 这 说 明 端 点 延 拓 方 法 能 够 有 效 地 抑 


制 最 小 二 乘 拟 合唱 


8 现 的 端 部 畸变 现象 。 
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图 1 UTI-UTC 观测 序列 及 其 UTIR-TAI 序 列 


Fig.1 Original UT1-UTC time-series and UT1R-TAI time-series 
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图 2 UTIR-TAI 序 列 及 其 端 部 数据 延 拓 点 


Fig.2 Original and extended UT1R-TAI time-series 
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3 UTIR-TAI 的 最 小 二 乘 与 基于 端 部 效应 改善 的 最 小 二 乘 拟 合 序列 


Fig.3 LS and ECLS fitting time-series of UT1R-TAI 
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4 UTIR-TAI 的 最 小 二 乘 与 基于 端 部 效应 改善 的 最 小 二 乘 拟 合 的 首尾 两 端 残 差 序 列 


Fig.4 LS and ECLS fitting residuals of UT1R-TAI at boundaries 


为 了 检验 端 部 效应 改善 的 最 小 二 


LS+NN 模型 和 ECLS+NN 模型 


Absolute Error，MAE) 作为 预测 性 能 


其 中 ， 已 ，O) 分别 为 了 点 的 UT1-UTC 预测 值 及 观测 值 ; 


本 文 总 共 作 了 300 期 的 预测 ， 
图 5 给 出 了 LS+NN 模型 和 ECLS+NN 的 UTI-UTC 预测 平均 绝对 误差 对 比 图 ， 表 1 给 
+ 在 不 同 预报 长 度 下 的 平均 绝对 误差 统计 结果 。 从 
WITHER, ECLS+NN 模型 的 预报 精度 相对 于 常规 LS+NN 模 


出 了 LS+NN 模型 和 ECLS+NN 模型 


表 ]1 可 以 发 现 ， 对 于 1 一 30d WER 


型 的 预报 精度 并 无 改善 


LS+NN 模型 ， 


MAE, = 


, (AM 30d 开始 ， 
预报 精度 最 大 提 


乘 外 推 模型 
I 对 UT1-UTC 作 1 一 360d 时 长 预报 ， 采 


评价 指标 ， 其 计算 公式 可 以 写 为 


| 


BI M =300 。 


外 推 模型 相 比 较 而 言 ， 端 部 效应 改善 的 最 小 二 
明显 的 优势 ， 同 时 也 从 侧面 反映 出 最 小 


影响 更 大 。 


(a)1~130d 


平均 绝对 误差 /ms 
a 


=== LS+NN 
---- ECLS+NN 


I 对 UT1-UTC 预报 的 改善 


, 


效果 ， 


分 别 利用 
平均 绝对 误差 (Mean 


(6) 


i 为 预报 长 度 M 为 预报 期 数 ， 


图 5 及 


ECLS+NN 模型 的 预报 精度 明显 优 于 
高 了 27%， 且 一 直 保 持 在 10% 以 上 ， 


平均 绝对 误差 /ms 
A 
o 


0 100 
预报 时 长 /d 


lá - 
' 


(b)130~360d 


200 250 300 
预报 时 长 /d 


图 5 LS+NN 模型 与 ECLS+NN 模型 的 UT1-UTC 预报 MAE 对 比 


350 


400 


Fig.5 MAE Comparison between the LS+NN and ECLS+NN models for UT1-UTC prediction 


K1 LS+NN 模型 与 ECLS+NN 模型 
Table 1 MAE statistics of UT1-UTC prediction obtained by the LS+NN and ECLS+NN models (unit: ms) 


预报 时 长 /d LS+NN 
1 0. 04 
10 0. 88 
20 2. 18 


ECLS+NN 
0. 04 
0. 87 
2. 16 


精度 改善 百分比 
0% 
1% 
1% 


I 的 UT1-UTC 预报 MAE 统计 (单位: ms) 


常规 


这 说 明 与 常规 最 小 二 乘 
乘 外 推 模型 对 于 UT1-UTC 中 长 期 预报 具有 更 
二 乘 拟 合 的 端 部 畏 变 现象 对 UT1-UTC 中 长 期 预报 的 


Chink ive" 


3 


30 3. 51 3. 44 2% 
60 7. 44 6. 49 15% 
90 12. 46 10. 59 18% 
180 26. 64 22.07 15% 
240 37. 39 27. 28 27% 
300 50. 88 38. 85 24% 


360 56. 78 50. 40 11% 
结论 


本 文 提 出 了 一 种 顾及 最 小 二 乘 拟 合 端 部 效应 的 UT1-UTC 预报 方法 ， 这 种 方法 与 基于 常 


规 最 小 二 乘 外 推 模型 的 UT1-UTC 预报 方法 区 别 在 于 ， 该 方法 在 进行 最 小 二 乘 拟 合 之 前 ， 首 
先 在 原始 序列 两 端 增加 统计 延 拓 数据 ， 然 后 再 对 数据 延 拓 后 的 新 序列 进行 最 小 二 乘 拟 合 ， 
其 目的 为 将 最 小 二 乘 拟 合 存在 的 端 部 畸变 搬移 到 新 序列 的 两 端 ， 从 而 抑制 最 小 二 乘 拟 合 序 
列 的 端 部 畸变 。 数 值 分 析 表 明 ， 最 小 二 乘 拟 合 的 端 部 效应 主要 影响 UT1-UTC 中 长 期 预报 结 


果 ， 


Me 


能 够 有 效 地 改善 最 小 二 乘 拟 合 出 现 的 端 部 畸变 现象 ， 与 常规 最 小 二 乘 外 推 模型 相 比 ， 基 于 
端 部 效应 改善 的 最 小 二 乘 外 推 模型 的 UT1-UTC 中 长 期 预报 精度 提升 明显 。 
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通过 在 观测 资料 的 两 端 增加 用 统计 学 方法 延 拓 的 外 推 数据 点 ， 然 后 再 进行 最 小 二 乘 拟 合 ， 


A least squares extrapolation model for UT1-UTC 
prediction method with consideration of the edge- 


effect 
Lei Yu'”, Zhao Danning'?*, Cai Hongbing!” 


(1. National Time Service Center, Chinese Academy of Sceinces, Xi an 710600, China; 2. Key Laboratory of 
Time and Frequency Primary Standards, Chinese Academy of Sciences, Xi’an 710600, China; 3. University 
of Chinese Academy of Sciences, Beijing 100049, China) 

Abstract: The prediction accuracy of UT1-UTC can be easily affected by the edge distortion 
of least squares (LS) fitting time-series, referred to as edge-effect in the data-processing domain, 
when periodic oscillations and residuals are separated by LS fitting. In order to alleviate the edge- 
effect, the original UT1-UTC time-series is first extended on both boundaries by using a time- 
series analysis model in this paper. A LS extrapolation model is then set up using the extended 
time-series. Finally UT1-UTC predictions are obtained by employing the combination of the edge- 
effect correlated least squares (ECLS) model and a stochastic predication technology such as 
neural network (NN). The numerical experiments demonstrate that the edge-effect can be 
noticeably alleviated with the developed method. In addition, the accuracy of the UT1-UTC short- 
term predictions is comparable with that by the conventional LS extrapolation-based prediction 


algorithm. However, the medium- and long-term predictions are significantly more accurate than 


those obtained by the proposed ECLS extrapolation-based prediction solution. 
Key words: UT1-UTC prediction; Prediction model; Least squares (LS); Edge-effect; Edge 


extension 


